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Section I 


INTRODUCTION 


A. Contract Objectives and Scope 


The objectives of this contract (i.e., Task Assignment 17) are to advance the hybrid 
phase stabilization (HPS) design procedures which were developed under Task 
Assignment 10 entitled, "Aeroservoelastic Stabilization Techniques for Hypersonic Flight 
Vehicles" for NAS 1-1 8763 (Ref. 1); to further develop residual response metrics as a 
basis for alternative structural stability specifications; and to develop strategies for 
validating HPS design and specification concepts in manned simulation studies. 

In this report, Task 10 and Task 17 are referred to as Phase I and Phase II, 
respectively. This report presents the results, conclusions and recommendations of the 
second phase of the contract. The Phase II work is a logical continuation of the Phase I 
work, both in development of improved flight control system designs and in refinements to 
the methods for assessing these designs. 

B . Technical Approach 

The motivation for the overall project began with a recognition that hypersonic 
vehicle (HSV), such as the National Aerospace Plane (NASP), can be expected to be 
particularly challenging with respect to the interaction of structural dynamics and flight 
control. Basically, this stems from the relatively low frequency structural modes inherent 
in the long, slender HSV configurations. The HSV model used to date has its first 
structural mode at about 2 Hertz (Hz). The low frequency modes create problems for 
flight control system (FCS) design which are compounded by the increased bandwidth 
requirements for relaxed static stability configurations such as the project model. 

Conventional treatment of these problems, i.e., gain stabilizatitMi in the FCS design 
using notch filters applied to selected structural modes, can, of course, be applied to meet 
existing gain and phase margin requirements of MIL-F-87242 (Ref. 2). However, as 
confirmed in the Phase I analysis, such notch filters applied to the low frequency structural 
modes of HSV inevitably create very high levels of effective time delay. This creates a 
significant problem in connection with manual flight control. The Phase I work indicated 
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the inherent difficulty in meeting the MIL-STD-1797 (Ref. 3) Level 1 effective time delay 
requirements using conventional gain stabili/ed FCS architectures. 

Reduction in the effective time delay to aid manual control motivated the primary 
flight control concept being pursued in this project - hybrid phase stabilization (HPS). 
HPS is simply the combination of conventional gain stabilization with phase stabilization. 
For this application, the essence of phase stabilization is manipulation of open loop zeros 
with respect to their companion poles so that there is a favorable increase in phase near each 
structural dipole. This creates a situation where closed loop stability will be desensitized 
with respect to feedback loop gain variations. As illustrated in Figure 1, this basically 
requires placing the zero below the pole. In the pitch rate-based FCS designs used in this 
work, zero placement is mechani^d by appropriate mixing of signals from a forward 
mounted and an aft mounted pitch rate gyro. The first task of the Phase II work was 
devoted to improving the Phase I HPS design. Significant improvements were achieved 
and these are discussed in Section II-B. 

Application of an unconventional FCS concept such as HPS requires careful 
consideration of the existing design requirements - not only to determine if they are 
adequate to insure safe designs, but also to insure that they do not unreasonably restrict 
application of new concepts. The second point is critical for HSV where unprecedented 
design challenges do not allow designers the luxury of overly conservative design 
requirements. The requirements must be as good and as well-founded as the FCS designs 
themselves! Thus, regardless of whether or not the HPS concept is ultimately the best 
design approach for HSV FCS, the intrinsic difficulty of this design problem requires 
careful examination of the design requirements. Such an examination has been part of this 
project from Phase I. 

Among the most fundamental design requirements are the gain and phase margin 
requirements of MIL-F-87242 (which have been carried forward from its predecessors 
MIL-F-9490D and earlier specs). These specifications are essentially traditional robustness 
requirements that have not been re-validated for the new HSV class. Furthermore, 
robustness specifications such as these are not explicitly linked to physical phenomena 
which can be re-examined in the context of a new aircraft type. Thus, as part of the 
Phase-II study, two phenomena of particular significance to the requirements for applying 
HPS to HSV have been examined. One, the impact of HPS on structural fatigue, is 
"aircraft-centered"; that is, it is independent of human pilots and passengers. The other, 
biodynamic interaction between the pilot, the FCS and the structural dynamics, is human- 
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centered (i.e., it would not occur in an automatic FCS mode). These analyses are 
discussed in Section IV. 

As noted above, the existing relevant design specifications, the gain and phase 
margin specs of MlL-F-87242, address robustness. However, application of HPS raises 
some new robustness questions, specifically the impact of uncertainty in structural mode 
zeros which are not explicitly addressed by the MIL-F-87242 specs. Thus in Phase n, the 
relative sensitivity of conventional and HPS designs to uncertainties in the structural model 
was examined. The results are discussed in Section III. 

The linearized dynamic model developed in Phase I for a power-on Mach 6 flight 
condition was used in Phase 11 with some extensions. Specifically, the original structural 
dynamic model was extended to generate a consistent representation of the normal 
acceleration (n 2 ^. This was done to replace the rough estimate of n^ used in Phase I which 
was not derived directly from the structural model. The improved nz model was used to 
revise the Phase I residual response metric (Section IV-A) and the analyses of aircraft- 
centered and pilot-centered phenomena impacted by HPS (Section IV). 

In Section II, the refinement of the nz model is discussed first, followed by the 
refinement of the Phase I HPS design. The remainder of the report is devoted to 
assessment of the final Phase II HPS design (referred to here as the "Phase Ila HPS" 
design). It will be seen that significant improvements were achieved in developing the 
Phase 11a HPS design, however, the assessments have given a clearer understanding of 
issues which should be further examined as well ideas for new approaches. These are 
discussed in Section V. Finally, the conclusions and recommendations for future research 
arc presented in Section VI. 
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Section II 


REFINEMENTS OF FLEXIBLE HSV MODEL AND HPS DESIGNS 


A. Refinement of Flexible HSV Model 

The flexible HSV model has been extended to provide normal acceleration per wing 
rotation at the pilot's station, fuselage station (FS) 84. The frequency response is shown in 
Figure 2. The development of the normal acceleration model is similar to that of the pitch 
rate model (Ref. 1). 

As shown previously in Phase I work, the HSV is modeled with uncoupled rigid 
body and flexible dynamics. These dynamics are combined as depicted in Figure 3. The 
derivation of the accelerometer output equation is presented in detail in Appendix C. 

B. Refinement of HPS Designs 

1 . Overview of Phase I and Phase II Designs 

At the conclusion of the Phase II work, four basic flight control system designs 
have been developed. These are designated as: 

• Gain Stabilized Reference Design 

• Phase I HPS Design 

• Initial Phase II HPS Design 

• Final Phase Ha HPS Design 

The final Phase Ila design is considered the best HPS design to date. The gain 
stabilized reference design was developed in Phase I and has been used through Phase II 
as the basic comparison with standard FCS design practice. A variant (referred to as "gain 
stabilized H") was considered in connection with the robustness analysis and is discussed 
only in Section III. Finally, reference is occasionally made to the "baseline" system. This 
system is not considered an actual flight control design because it consists of only the 
superaugmented pitch loop designed in Phase I based on the rigid body dynamics without 
treatment of the structural dynamic modes. 
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Figure 3. Aircraft Dynamics Modeling 


All four designs, along with the baseline system, can be represented by one generic 
block diagram as shown in Figure 4. The Gc element is the basic rigid body control law 
which is the same in all designs and would constitute the final design if there were no 
structural modes. The gain stabilized reference design uses only the signal from the 
forward pitch rate gyro (FS 84) passed through a notch filter (Gsi, placed near the first 
structural mode) and a higher frequency lag (Gs 2 ). The frequency response of the gain 
stabilized design is shown for reference in Figure 5. All three HPS designs are based on 
the blending of signals from forward and aft pitch rate gyros suitably filtered to position the 
zeros associated with the lower frequency structural dipoles. The differences in the three 
HPS designs are in the Gsi and Gs2 elements. The evolution of these elements are 
discussed in the following subsections. The block containing the gain function 
l/(Ka + Kf ) was included to insure that the gyro signal blending would not change the 
rigid body crossover frequency and distort the closed loop short period mode. 

2 . Review of the Phase I HPS Design 

As noted previously, the HPS concept is to combine conventional gain stabilization 
with phase stabilization in a way that exploits the advantages of both schemes and 
accommodates their disadvantages to produce a superior design. However, a wide variety 
of design requirements must be considered in the design tradeoffs. In reality, and of 
particular concern for HSV designs, these requirements must include much more than the 
gain and phase requirements of MIL-F-87242. These additional "implied requirements 
have been an integral part of this study from the beginning and are discussed in more detail 
in Section IV. 
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Figure 4. Flight Control Systran Block Diagram 
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The requirement motivating HPS is the need to reduce effective time delay 
associated with the notch and lag filters used in conventional gain stabilized designs. As 
will be seen, HPS provides significant time delay reductions as well as improvements in 
stability. The disadvantages of HPS are associated with increased "residual 
(uncommanded) response" arising from the structural modes. In addition, because phase 
stabilization depends on the proper placement of open loop structural zeros, there is 
concern about the robustness of phase stabilized designs with respect to uncertainty in these 
zeros. However, conventional designs have their own robustness problems if notch filters 
are not placed and scheduled properly. The robustness issue is discussed in Section HI. 

The initial investigations of Phase I led to the approach of phase stabilizing only the 
first structural mode and using a lag filter (but no notch filter) to gain stabilize the higher 
frequency mode. This was truly a hybrid architecture where application of phase 
stabilization only to the first structural mode maximized the reduction of effective time delay 
(by eliminating the low frequency notch). It also minimized robustness problems on the 
assumption that at least the first structural mode zero could be defined with reasonable 
certainty. 

The potential of this strategy can be seen in Figure 6, a special root locus plot in the 
region of the first structural mode. Here, the basic (rigid body) pitch loop is closed (with 
only the Gc equalization) using the pitch rate signal from one of four different gyro 
locations for each closure. The range of zero locations seen in Figure 6 suggests that, by 
blending signals from several gyros, an effective zero could be located to insure phase 
stability of the first structural mode. This was done in the Phase I HPS design as 
documented in Ref. 1. The Gf^ and Gfj filters in the forward and aft pitch rate signal 
paths, as shown in Figure 4, allow the effective zero to be placed at will. Figure 7 shows 
a system survey of the pitch loop closure for the Phase I HPS design. 

The first level of assessments for all of the designs are related to stability. These 
are summarized in Table 1 in terms of gain, phase and time delay margins. The last 
column is effective time delay which is not strictly a stability issue for the flight control 
system loop closure, but is of primary concern for the human pilot's closed loop manual 
control. The cells in Table 1 are hatched where the corresponding stability metrics do not 
meet the MIL-F-87242 requirements. (Note there is no MIL-spec requirement on delay 
margin). It can be seen that all of the designs fail to meet the low frequency gain margin 
specification. This occurred because no specific effort was made to meet this specification 
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Figure 6. Pitch Loop Qosures at the Ffrst Bending Mode for Four Sensor Positions 










Table 1. Coaq>arison of Stability Metrics 
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which is not of particular concern in this study. This issue is discussed briefly in 
Appendix A. 

The major problem of the gain stabilized reference design is the expected large 
effective time delay which grcady exceeds the MIL-F-87242 maximum. The initial Phase I 
design resulted in a significant reduction in the effective time delay below the MIL-spec 
limit. Both the gain stabilized design and the Phase I HPS design also fail to meet the low 
frequency phase margin requirements, although the Phase I HPS design comes 
significantly closer. The stability assessment of the Phase II HPS designs will be 
discussed in the following subsections. 

Stability can also be viewed by comparing the closed loop damping ratios of the 
structural modes. This only provides some insight into the relative effect of the different 
designs on structural damping. Phase stabilization should generally increase closed loop 
damping and this is seen, for the first structural mode, when the Phase I HPS design is 
compared to the gain stabilized reference design. However, as shown in Table 2, the 
damping ratios for the second and third structural modes are actually reduced by the 
Phase I HPS design. This is consistent with the fact that no attempt was made to extend 
phase stabilization to higher modes to enhance the damping. Assessing the overall 
significance of these various damping ratios requires measures of residual response which 
will be addressed in Section IV. In terms of residual response, however, the Phase I 
HPS design was judged significantly deficient (Ref. 1) and this largely motivated the 
Phase II HPS refinements discussed next. 


Table 2. Comparison of Flexible Mode Damping Ratio 


Design Case 

Flexible Mode Qosed Loop Damping Ratios 

1st Mode 

2nd Mode 

3rd Mode 

4th Mode 

5th Mode 

Gain Stabilized Ref. 

0.00097 

0.01940 

0.01070 

0.00942 

0.00998 

Phase I HPS 

0.00471 

0.00629 


0.01020 

0.09700 

Initial Phase II HPS 

0.00781 

0.01640 

0.00808 

0.01080 

0.01040 

Final Phase Ila HPS | 

0.00842 1 

0.02870 1 

0.01490 

0.01040 

0.01310 
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The Initial Phase II HPS Design 


The essence of the Phase I HPS design process was to position the zero more or 
less directly below its companion pole in the first flexible mode dipole. The separation of 
the pole and zero is set by adjusting the lead in the filter Gf, in the aft gyro signal path. In 
the Phase I design this lead was set at 12 rad/sec. Among the first refinements considered 
in Phase II was adjustment of the first flexible dipole separation using the Gf, filter. A 
survey of the sensitivity to the zero of this flexible dipole, beyond what was done in 
Phase I, revealed that there was very little potential for design improvement at this level 
(Ref. 4). 

Thus, the Phase I HPS design was examined at a more fundamental level in a 
search for improvements. The Bode root locus of Figure 7 indicates that the second 
structural mode is critical in the Phase I HPS design. (The dots in the Bode root locus are 
open where the closed loop roots are unstable indicating a gain margin of about 8 dB at the 
second mode). However, examination of the overall Figure 7 system survey reveals some 
basic topological anomalies in the root loci of this design. These are seen most readily in 
the conventional root locus of Figure 7 where the locus departs directly to the right half 
plane from the second flexible mode. The companion zero is serviced by a locus emerging 
from the real axis. Thus, in a sense, this design causes the rigid body dynamics to couple 
to the structural dynamics in a way that is not particularly desirable. Furthermore, 
examination of the third structural dipole shows that, while it is phase stable in 
configuration (i.e., the zero is below the pole), the local root locus does not reflect this. It 
is not concave into the left half plane based on the Figure 1 ideal. 

This insight focused development on the potential for phase stabilizing the next 
several structural naodes as well as the first. Figure 8 displays the root loci near the second 
and third flexible modes. As in Figure 6, the basic (rigid body) pitch loop is closed (with 
only the Gc equalization) using the pitch rate signal from just one of four different gyro 
locations for each of four closures. It can be seen that, as with the first flexible mode, 
blending the FS 84 and FS 1050 pitch rate signals has some apparent potential for phase 
stabilizing the second and third flexible mode. 

Thus, the basic pitch blending used in the Phase I HPS design was examined at the 
second and third flexible modes as was done for the first mode in the Phase I HPS design 
(Ref. 1). As shown in Figure 9, the blending of the FS 84 and FS 1050 sensor signals, 
using the established blending filters but with no lag in Gsi, was treated as a root locus 
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a) Second Structural Mode 



b) Third Structural Mode 


Figure 8, Root l^ocus of 2nd and 3rd Bending Modes for Each of Four Sensor Positions 
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problem. It can be seen that phase stable dipoles can be created using the same Ka and Kf 
gain of the original Phase I HPS design. This established the fact that the topological 
anomaly seen in the Figure 7 root locus was created by the lag in Gsi. 

Having established that the next several flexible modes could also be phase 
stabilized, the lag filter in the element Gsi was replaced with a lead filter to achieve 
extended phase stabilization. To set the frequency of the lead, a survey was made in which 
the lead frequency was varied from 20 to 40 rad/sec. A lead frequency of 40 rad/sec was 
selected for two reasons. First, it provided the maximum potential damping of the four 
phase stabilized modes. Second, for the phase unstable fourth flexible mode, the selected 
lead filter provided the maximum gain margin. 

Figures 10 and 1 1 show the system survey plots. From the Bode root locus plot in 
Figure 10, it can be seen that the 8 dB gain margin requirements have been exceeded by at 
least a factor of two. Note that all modes, except the fourth flexible mode, are phase stable. 
The Bode plots of Figure 1 1 show the gain and phase margins tabulated in Table 1. 

This design constitutes the initial Phase II HPS design which is compared to the 
two Phase I designs in terms of the stability metrics in Table 1. The initial Phase II HPS 
design shows significant reduction in the effective time delay, an important result since this 
motivates the HPS concept. This design does not have the low frequency phase margin 
problem that the earlier designs exhibited. However, it does fail to meet the MIL-spec 
phase margin at the third flexible mode. This mode is phase stable so gain uncertainty is 
not a problem, but phase uncertainty (such as from unmodeled higher frequency modes) 
remains a concern. 

The closed loop damping ratios of the structural modes with the initial Phase II 
HPS design are compared to those in the two earlier designs in Table 2. It can be seen that 
the damping ratios of the initial Phase II design are generally somewhat higher (first and 
second modes) or about the same as the HPS Phase I. However as noted above in 
connection with Table 2, the overall impact of these damping changes cannot be assessed 
without addressing residual response metrics as will be done in Section IV-A-3. But as 
shown in Figure 12, the closed loop pitch rate response (at the pilot's station, FS 84) to a 
unit step input does give some indication of the residual response. Comparisons of 
Figures 12 and 13 indicate that the residual response of this HPS design is similar to that 
of the HPS Phase I design, but both of these designs have more residual response than the 
gain stabilized reference. The dominant 2 Hz residual oscillation seen in the HPS designs 
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Figure 12. Closed Loop Pitch Rate Response at the Pilot's Station to a Unit Step Input 

for the Initial Phase II HPS Design 
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Figure 13. Response Comparison of Pitch Rate to Unit Pitch Rate Command 
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arises from the first structural mode which is, of course, notch filtered in the gain stabilized 
design. 


In general, phase stabilization does have the potential to significantly reduce 
residual response by increasing the damping (as distinct from reducing magnitude of the 
response with a notch filter). However, this potential is not realized in the first two HPS 
designs. Insight into this can be obtained from examination of the root loci in Figures 10 
and 14 for the initial Phase II HPS design. It can be seen that, for each of the structural 
modes, the actual damping ratio achieved is considerably less than what could be achieved 
if the loop gain was increased. It is important to note that a simple increase in the loop gain 
cannot be made to improve the flexible mode damping without seriously degrading the 
closed loop short period (rigid body) dynamics. Design alternatives to resolve this conflict 
lead to the final Phase Ila HPS design discussed in the next section. 

4. The Final Phase Ila HPS Design 

The basis of the final Phase Ila HPS design was to boost the loop gain in the 
flexible mode frequency range in order to increase the flexible mode damping without 
changing the loop gain in the region of the rigid body crossover which sets the closed loop 
short period dynamics. This suggested inserting a filter into the pitch loop, with unity gain 
at low frequencies and higher gain in the structural mode region. The first order lead-lag 
filter shown in Figure 4 was used in the Gs2 element to accomplish this. An alternative of 
using a second order lead-lag filter was also examined (Ref. 4) but the first order filter 
selection proved superior and this Phase Ha HPS design is the only configuration that will 
be discussed here. 

The system survey for the Phase Ila HPS design is shown in Figure 15. The zero 
of the lead-lag filter was placed at 10 rad/sec just below the first flexible mode frequency. 
This flattens the Bode asymptote which provides increased gain and thus increased 
damping to the phase stabilized flexible modes. The pole was placed at 25 rad/sec to gain 
stabilize the higher modes consistent with the HPS concept. This lag break frequency was 
selected to provide the required 8 dB gain margin at the fourth flexible mode - the only 
mode that is not phase stable. 

From Table 1 , it can be seen that the final Phase Ila HPS design results in the 
lowest effective time delay (36 msec) of the four designs. This is notably lower than the 
baseline system which has no added structural filters. This, of course, results from the 
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Figure 15. System Survey of Phase Ila HPS Pitch Loop Closure 
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lead introduced to increase the structural mode damping. It is worth noting that there are 
alternatives to phase stabilization which would also effectively introduce lead. For 
example, various pitch acceleration estimation schemes exist that are based on pitch rate 
measured at a single fuselage station. This suggests a more complete investigation of the 
relationship between these two concepts as part of future design refinements. 

The gain and phase margins of the final Phase Ila HPS design are shown in 
Figure 16. There is a phase margin subtlety which applies to other designs and certainly 
should be noted for the final Phase Ha HPS design. Near structural modes, for example 
the second mode, the magnitude peaks cross the 0 dB line - an unconventional situation. 
Thus, there are phase margins defined at each of the zero crossings near the second flexible 
mode peak (-166.7® at 36.3 rad/sec and 79,8® at 37.8 rad/sec). The negative phase 
margin deserves some comment. The negative phase margin can be converted to a positive 
value by adding 360®. Hence, the phase margin -166.7® is equivalently 199.3®. By 
conventional reasoning, the smaller positive margin (i.e., 79.8®) is the more critical. In 
constructing Table 1, only the smallest positive phase margin of a modal pair is listed. 
However, for HPS designs, the customary thinking needs to be examined. This is done in 
Appendix B where it is noted that large positive phase margins do not necessarily imply 
robustness. 

Table 1 indicates that the final Phase Ila HPS design meets all applicable MIL-spec 
requirements except, as with the initial Phase II HPS design, the phase margin at the third 
flexible mode. This phase margin is actually lower for the final Phase Ila HPS design than 
the initial Phase II design. As noted above, this mode is phase stable which reduces 
concern for gain uncertainty but not phase uncertainty. 

The conventional root locus of Figure 15 shows that the damping of the first 
flexible mode has been slightly increased compared to the initial Phase II design 
(cf. Figure 10). Table 2 shows similar increases for the other phase stabilized flexible 
modes and, in fact, the final Phase Ila HPS design exhibits the highest structural damping 
ratios except for the fourth flexible mode where the initial Phase II HPS design produces 
higher damping. As noted previously, in connection with Table 2, no overall conclusions 
can be drawn from these damping ratios, llie closed loop pitch rate response to a unit step 
pitch rate command shown in Figure 17 suggests some residual response improvement 
with respect to the initial Phase II HPS design, but the gain stabilized design still appears 
to have the advantage in residual response over all HPS designs. The significance of these 
differences will be addressed in Section IV. 


26 




Vigfut 16. Open Lcx>p Frequency Response of Rnal Phase Ila HPS Design 
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Figure 17. Closed Loop Pitch Rate Response at the Pilot's Station to a Unit Step Input 

for the Phase Ila HPS Design 
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Finally, it should be noted that the evolution of the HPS design has led to a final 
Phase Ila HPS design which is clearly unconventional. Specifically, the lead-lag filter in 
Gs 2 produces an effect directly opposite to the notch and lag filters used in conventional 
gain stabilization (which can be viewed as an "inverse notch" filter). The premise here is 
that the improved flexible mode damping will reduce the residual response to acceptable 
levels. This is certainly consistent with the original idea that HPS design is a matter of 
finding the best mix of gain and phase stabilization. This in turn suggests a formal optimal 
control approach. However, as will be further supported by the developments in 
Section IV, the complex and specialized design requirements of this problem require more 
attention at this point than the mechanics of optimization. 

The very low value of effective time delay achieved in the final Phase Ha HPS 
design (36 msec) implies a potential for further design optimization. Because there is no 
benefit for grossly exceeding the MIL-spec requirement, the effective time delay could be 
relaxed to tune the design for other requirements. Again, it is the specialized requirements 
forHSV design that are most critical at this stage. 
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Section III 


SENSITIVITY OF HPS DESIGNS TO STRUCTURAL UNCERTAINTIES 


Variations in structural mode parameters can have a significant destabilizing effect 
on a control system that would otherwise be considered a good design with proper and 
adequate compensation of structural mode interaction. This section examines the sensitivity 
of HPS designs to structural uncertainties. Results of two gain stabilized designs aie also 
included for reference and comparison. 


A . Modeling 


The flexible vehicle model was developed using FAMUSS, an aeroservoelastic 
modeling tool widely used at MCAIR (Ref. 5). The FAMUSS model can either be 
represented by a transfer function (pole- zero) form or the standard state space form. For a 
single- input single-output system, the generic state and output equations, respectively, are 


i = As, + B u 
y = Cs + Du 


where 


r Ai 


A = 


A2 


0 



B = [ 0 I • 


( 1 ) 


C = [ Cl C2 • • Cj ] 


D = d . 

For the analysis that follows, three structural mode parameters are of primary 
concern: natural frequency (o)s)> damping ratio (Cs)t ^ud nKxie shape (<^s)- Th® FAMUSS 
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model lends itself nicely to incorporation of structural uncmainties. This becomes apparent 
by observing FAMUSS realization of A and C. Specifically, natural frequency and 
damping ratio are uniquely defined in A while mode shape is contained in C. Using the 
subscript "i" to indicate the ith structural mode, 

1 

2 Cs,i ®s,i 

( 2 ) 

This realization allows direct substitution of the nominal values of natural 
frequency, damping ratio and mode shape with some new values that represent variations 
of these structural parameters. The new or corrupted values can be expressed as follow: 

“..I = (l±Ao)l»..i 

c.,i = (>± 

^ (3) 

where Aa>, A^, and A^ are the normalized uncertainty levels for natural frequency, damping 
ratio, and mode shape, respectively. Note that these parameters can be varied 
independently. 

Figure 18 shows the frequency response characteristics of the flexible vehicle 
model alone. The mode shape uncertainty changes the slope while natural frequency 
uncertainty primarily shifts the structural mode up and down. Although not obvious from 
the figure, damping ratio uncertainty changes the peaks and valleys of the frequency 
response. 

The pole-zero characteristics of the flexible vehicle model can be seen by expressing 
the input-output relationship in terms of the system matrices, i.e., 

y CAdj(sI- A)B + det(sI-A)D ... 

u " det(sI-A) • 


Ai = 


0 


-w 


S,l 


U i — <|>s,i 
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MODE SHAPE VARIATION 



FREQUENCY VARIATION 



DAMPING VARIATION 



Figure 18. Frequency Responses of a Flexible Vehicle Subject to Structural Uncertainties 
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Equation 4 suggests that while mode shape uncertainty affects only the zero 
location of a structural mode, natural frequency and damping ratio uncertainties alter both 
pole and zero locations of a structural mode. 

The uncertainty levels associated with the structural dynamics of the Space Shuttle 
Orbiter (Ref. 6) arc listed in Table 3. These uncertainty levels were chosen for this 
contract since the Space Shuttle Orbiter is the only operational hypersonic vehicle to date. 
The same set of data was also used in the Modem Aerospace Vehicle Robust H-Infinity 
Control (MAVRIC) program, published in Ref. 7. 

B. Analysis 

The design evaluation space includes two gain stabilized and two HPS designs. The 
Erst gain stabilized design, labeled GS-I, is the gain stabilized reference design shown in 
Figure 4. The second gain stabilized design, labeled GS-II, employs a more conservative 
approach to gain stabilization than GS-I. The two HPS designs arc Phase I and final 
Phase Ila HPS designs. Specific characteristics identifying each design arc summarized in 
Table 4. It must be emphasized that each design will be evaluated against its own 
rcfercnce case, i.e., one that has no variation in any of the three stractural mode parameters 
of interest. 

Since the damping ratio uncertainty level was rather small, the uncertainty level was 
arbitrarily increased from ±2% to ±15%. The effect of this damping ratio uncertainty alone 
was determined. It was found that both HPS and gain stabilization were insensitive to 
variations in damping ratio. Thus, damping ratio was held at nominal value (i.e., = 0) 

in subsequent analyses. 


Table 3. Structural Mode Uncertainty Levels 


PARAMETER 

UNCERTAINTY LEVEL (%) 

Frequency 

±10 

Mode shape 

±15 

Damping ratio 

±2 
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Table 4. Design Evaluation Space 


GS-I 

GS-II 

PHASE I HPS 

PHASE HA HPS 

U.se of ?wd sensor 
signal only 

Use of^ fwil sensor 
signal only 

Blending of fwd and 
aft sensor signals to 
phase stabilize first 
mode 


Notch and low pass 
filters to gain 
stabilize structural 
modes 

Notch and low pass 
filters to gain 
stabilize structural 
modes 

First order lag to 
gain stabilize higher 
frequency structural 
mo^s 

First order lead to 
enhance loci of 
higher structural 
modes 

rSBM 

Aimed to achieve 
8 dB peak clearance 

— 

First order lead-lag 
to increase structu^ 
damping 


The stability margin results of the two gain stabilized designs are shown in 
Tables 5 and 6. As expected, the low frequency gain and phase margins as well as the 
high frequency gain margin (at whfi) compare well with the reference case in the control 
bandwidth region (with frequency < 10 r/s). 

In the structural mode region (with frequency > 10 r/s), as the structural modes 
move up and down due to variation in frequency, the notch filter frequency no longer 
coincides with the structural mode frequency, causing insufficient gain attenuation at the 
first mode. This in turn causes the magnitude of the first mode to cross over the zero dB 
line resulting in an additional phase margin. See Figures 19 and 20. Note from Table 5 
or Figure 19 that another phase margin exists at the second mode for GS-I. This can be 
explained by recognizing that instead of meeting the more conservative 8 dB peak 
clearance, the GS-I design goal is to meet only the 8 dB gain margin definition. Table 5 
indicates that several gain margins at o)hf 3 (in the vicinity of the first mode) are well below 
the 8 dB spec. In addition, negative phase margins exist at the first mode for the cases with 
+10% frequency variation. Unlike Phase Ila HPS design (cf. Section II-B-4), the 
negative phase margin here implies instability. A detailed discussion of interpreting 
negative phase margin is provided in Appendix B. 

The sensitivity results of GS-Il is shown in Table 6. When compared with GS-I, 
GS-II appears to be more robust in meeting the gain margin requirement. This happens at 
the expense of a reduced low frequency phase margin, however. The phase margins at the 
first mode for GS-II do not quite meet the 60” spec. As expected, the results suggest that 
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Figure 19. Frequency Reqxmse of GS-I Design with -10% Frequency VariaticHi 
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Figure 20. Frequency Response of GS-H Design with +10% Frequency Variation 




gain stabilization, in general, is sensitive to variations in frequency, particularly at the first 
mode. 


Tables 7 and 8 show the stability results of Phase I HPS and Phase Ila HPS, 
respectively. With the exception of a few cases where the gain margins actually fall below 
the desired value of 8 dB, the results do not show any significant cause for concern. 
Generally, HPS exhibits robust characteristics at the first elastic mo(k. This is illustrated in 
Figure 21. Conceptually, the objective of phase stabilization is to position the zero of a 
structural pole- zero pair below the pole. As shown in the figure, even in the presence of 
structural uncertainties, the right half plane zero can generally be placed somewhere below 
the pole. For higher frequency, the uncertainty regions of the pole and zero can potentially 
overlap. This obviously is a concern and thus a primary reason of the use of HPS instead 
of phase stabilization alone. Finally, like conventional gain stabilization, HPS is sensitive 
to variations in frequency, but it generally occurs at a much higher frequency. 

Next, the root loci of the four designs are compared to its own reference case. It 
was observed that, when the frequency was varied, the root loci of GS-I and Phase I HPS 
were not maintained. See Figures 22 and 23. This would probably affect the loop gain 
sensitivity and ultimately the stability margins. 

At this point, it would be unfair to claim that HPS is more robust than gain 
stabilization or vice versa. The point to make here is that each stabilization technique, 
whether it is gain stabilization or HPS, has its own unique problems and robusmess 
characteristics. Tradeoff exists even within each individual class of stabilization 
techniques. The bottom line is that robustness or lack of it would likely not be the primary 
reason for prohibiting the use of hybrid phase stabilization. 
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Table?. Sensitivity of Phase I HPS Design to Si 
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Table 8. Sensitivity of Phase Ha HPS Design to Structural Ui 


Gain Margins (dB) 

U. 
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3 
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7.24 

9.45 

9.29 

8.13 

10.26 

7.70 

6.47 

8.76 

Low 

Freq 

1 
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-5.83 

-5.81 

-5.81 

CM 

00 

iri 

1 
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00 

iri 

1 
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00 

tri 
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Phase Margins (Deg) 
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CO 

CNJ 

27.01 

24.79 

29.94 
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CO 

CD 
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Figure 21. Conceptualization of Phase S^bilization in the Presence of Uncertainties 


43 




rrr* 



^se ♦10% Frequency Vari 




• • • • ^ 

o o o O o 

00 ^ 


AyVNlOVl^l 


45 


Figure 23. PhaselHPS Root Locus Compaii 
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Section TV 


REQUIREMENTS DEVELOPMENT AND SYSTEM ASSESSMENT 


A. Residual Response Metrics and Assessments 

Based only on qualitative comparisons of pitch rate responses to step pitch rate 
commands (Section II), all of the HPS designs showed evidence of residual oscillations 
associated with the first structural mode which exceed those of the gain stabilized reference 
design. This is not surprising since the HPS designs do not have the notch filter at the first 
flexible mode used in the gain stabilized design. If this were the only issue, there would be 
no reason to consider the HPS designs further. There are of course other issues. Primary 
issue anoong these is the unacceptable effective tirne delay of the gain stabilized design, the 
problem which motivates HPS. As noted previously, designers will not have the luxury of 
overly conservative flight control design for HSV. Therefore, the tradeoffs between time 
delay, robustness and residual response will have to be made very carefully. All of the 
requirements that effect these tradeoffs will have to be well defined, traceable to physical 
principles and well validated for HSV. The greatest need is in the area of residual resptmse 
requirements which are nearly nonexistent for HSV. 

A residual response metric was developed in Phase I but it was recognized to be 
very questionable, in part because it was based on an inadequate estimate of the normal 
acceleration at the pilot's station. This problem has been remedied and an improved 
"Phase H" residual response metric will be reviewed later. However, even this Phase II 
metric is inadequate because it is still based on an ad hoc frequency dependent weighting. 
To address this limitation, studies of relevant aircraft-centered and human-centered 
phenomena were made to provide the foundation for more rational residual response 
metrics. These are presented in Sections IV-B and IV-C. 


1 . Normal Acceleration Frequency Response 


As discussed previously in Section II-A, the structural dynamics model was 
extended in Phase II to properly compute the normal acceleration at the pilot's station. 
Figure 24 shows the magnitude plots of the pilot's station normal acceleration to pitch rate 
command for: (a) the baseline system (no stmctural mode treatment); (b) the gain stabilized 
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Figure 24. Noirnal Acceleration Magnitude Plots at the Pilot's Station 
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reference design; and (c) the final Phase Ila HPS design. The corresponding transfer 
functions are given in Table 9 in factored form. As expected, the two structural mode 
treatments each show clear reductions in normal acceleration at the structural modes when 
compared to the baseline case. 

From Figure 24, the gain stabilized design shows somewhat lower acceleration 
peaks than the Phase Ha HPS design at the first, third and fifth structural modes but not at 
the second and fourth. Thus, a clear conclusion about structural response cannot be made 
from qualitative comparison of frequency response. The question ultimately requires that 
some frequency dependent weighting be applied to the response. Defining an appropriate 
weighting is, of course, the essence of the development of a residual response metric. 

2 . Normal Acceleration Step Response 

Figure 25 shows the normal acceleration at the pilot's station (in g s) in response to 
a 1 °/sec pitch rate command step comparing the gain stabilized reference and the final 
Phase Ila HPS design. Figure 26 presents the same comparisons over a shorter time 
interval. The generally slow (rigid body) response of acceleration to pitch rate is a 
consequence of the pure pitch rate command stability and control augmentation system 
(SCAS) design coupled with the unconventionally low values of I/T 02 typical of 
hypersonic flight (Ref. 8 and 9). This characteristic is quite important in hypersonic flying 
qualities and has an impact here as will be discussed later. 

The step response calculation has the benefit of applying an intrinsic weighting to 
the modal responses but only subjective assessments can be made by comparing the 
responses of the two designs. On this basis, however, no major differences are indicated 
between the Phase Ila HPS design and the gain stabilized reference design. The residual 
response of all designs appears dominated by the 53 rad/sec signature of the third 
structural mode which is consistent with the frequency responses shown in Figure 24. For 
the HPS designs (Figures 12, 13, and 17), the significance of the third structural mode in 
the O/p step responses contrasts with the 2 Hz residual oscillation from the first mode seen 
in the pitch rate responses. This arises because the higher frequency modes are accentuated 
by the inherent pitch rate to t'zp due to gyro offset from the c.g. contribution. 
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Table 9. Nonnal Acceleration Transfer Functions at the Pilot's Station 


Baseline 

4.98E+07(-7.40E-04)(.44)(2.31)[-.04, 2.7^[.12, 12.34] 
”Zp [-2.78E-03. 25.94][.38. 52.41 ][-e.69E-03, 67.99][.049, 74.66] 

^ “ (0)(. 16)[. 70. 1.46]{2.31)[4.03E-03, 12.50](25.99)[.02, 36.68] 

[2.21 E-03, 53.52][.012, 64.71 ][9.99E-03, 77.03][.51, 272.93] 


Gain Stabilized Reference Design 


7.96E+10(-7.40E-04)(.44)(2.31)[-.042. 2.79K.12, 12.34][.01, 12.5] 
^Zp [-2.78E-03, 25.94][.38. 52.41 ][-8.69E-03. 67.99][.049, 74.66] 

“ (0)(.16)[.82, 2.05](2.31)[.66, 1 1 .00][9.76E-04. 12.6^(18.60) 

[.019, 37.13][.53, 37.77][.011, 53.1^l9.42E-03. 64.71] 
(9.98E-03, 77.02](.51, 272.89] 


Final Phase Ha HPS 


4.98E+07(-7.40E-04)(.44)(2.31)[-.042. 2.79][.12. 12.34](25)3 
^2p [-2.78E-03. 2S.93][.38, 52.41 ][-8.69E-03, 67.99][.05, 74.66] 

^ ' (0)(.14)[.67, 1.20](2.31)[8.42E-03. 12.60](25)(25.30)[.96.31.20] 

(2.87E-02. 37.25][.015, 54.79](.010, 64.32H.013, 77.55J[.53, 253.23] 
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Figure 25. Responses of Normal Acceleratiwi at the Pilot’s Station to a Unit Pitch Rate 

Step Command - 5 Second Interval 
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Figure 26. Responses of Normal Acceleration at the Pilot's Station to a Unit Pitch Rate 

Step Command - 2 Second Interval 







3. 


Residual Response Metric Revisions 


The acceleration residual response metric developed in Phase I (Ref. 1) is 
quantitative and based on a physically motivated, but ultimately arbitrary, frequency- 
dependent weighting. In Phase I, a crude estimate of normal acceleration at the pilot's 
station had to be used but this was replaced with the much better acceleration transfer 
function now available (Section IV-A-1) to form an improved "Phase-II" metric. The 
other elements such as the weighting function and the random disturbance injected at the 
elevator surface remains unchanged from the Phase I development. The Phase I and H 
residual response metrics are summarized in Figure 27. 

The Phase I and Phase II residual (acceleration) response metrics are compared in 
Table 10. The values of the metric (rms acceleration) are generally an order of magnitude 
higher with the improved model which simply reflects the inadequacy of the Phase I 
normal acceleration model at high frequencies. Sec the frequency response comparison in 
Figure 27. However, the relative ranking of the designs does not change with the 
Phase II nnctric except that the Phase Ila HPS design appears equal to rather than superior 
to the gain stabilized design. 

As noted above, even with an adequate normal acceleration model, the Phase n 
residual response metric must be questioned because it is based on a frequency domain 
weighting function which is ultimately arbitrary and is not based on relevant physical 
phenomena. Thus, two important areas were studied to develop deeper insight into the 


Table 10. Comparison of Residual Response Metrics 


System 

Phase I Metric 

IHiase n Metric 

o„,(g/deg) 


On^(g/deg) 


"Baseline" 

0.0872 

1.000 

0.750 

1.000 

Gain Stabilized Ref. 

0.0587 

0.673 

0.350 

0.476 

Phase I HPS 

0.0609 

0.698 

0.488 

0.651 

Hnal Phase Ila HPS 

0.0441 

0.506 

0.350 

0.467 


53 



Disturbance 

Spectral Augmented Residual Response 

Filter Aircraft Weighting Filter 



54 


Figme27. Residual Response Metric 





significance of residual response and to search for more rational and physically traceable 
response metrics. These can be broadly partitioned inu> two groups: aircraft-centered issue 
which is independent of the human pilot or passengers, and human-centered issue which is 
dependent of the human pilot. The aircraft-centered issue of structural fatigue is discussed 
next in Section IV-B. The discussion of the human-centered problem of biodynamic 
vibration transmission to the stick follows in Section IV-C. 

B. Aircraft-centered Requirements 

Metal fatigue is a phenomenon experienced by structures when exposed to repetitive 
load cycles. These loads cause fatigue cracks which grow until fracture of the structural 
member occurs. The fatigue life of a structure is dictated by the number of cycles that 
cause fracture. Factors which affect fatigue include environmental conditions ( that is, 
elevated temperature exposure) and stress concentrations (e.g., fastener holes, cutouts, 
etc.). 


An assessment was performed to compare the fatigue life of a structure under loads 
caused by conventional gain stabilization and HPS. The time history of the normal 
acceleration at the pilot's station for a conventional gain stabilized design is shown in 
Figure 28. A similar time history plot for Phase Ila HPS is shown in Figure 29. These 
accelerations are assumed to be the loads that cause corresponding cyclic stress in the 
structure. 

The material selected for the structural fatigue life analysis or assessment is 
titanium 6-4. This is a conventional titanium alloy used extensively for elevated 
temperature conditions. The titanium 6-4 fatigue properties required in the analysis were 
obtained from MCAlR's in-house material database and the military handbook (Ref. 10). 
Structural failure is defined by the predicted initiation of fatigue cracks. The tool used for 
the analysis is SPECCI, a MCAIR fatigue crack initiation prediction software (Ref. 11). 
Two (low and high) temperature settings were studied: room temperature and 600°F. 

Results of the room temperature analysis for the gain stabilized design are shown in 
Figure 30. Stress on the vertical axis is in units of kilopounds per square inch (ksi) while 
the fatigue life on the horizontal axis is in number of maneuvers which cause structural 
failure. The maneuver for the gain stabilized design is the 15 second time history data as 
shown Figure 28. Similarly, the maneuver for the hybrid phase stabilized design is the 
15 second time history data shown in Figure 29. 
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Figure 29. Time History of Normal AcceleratiOT Response at the Pilot's Station to a Unit 
Pitch Rate Step Command for Rnal Hiase BLa HPS Design 



GAIN STABILIZATION FATIGUE AT BOON TENP. 

HYBRID PHASE STABILIZATION FATIGUE AT ROOM TEMP 
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Rgure 30. CompaiiscKi of Predicted Fatigue Life at Room Temperature 




As illustrated in Figure 30, gain stabilization results in longer fatigue life than 
hybrid phase stabilization. This difference can translate into structural weight savings. For 
example, if a vehicle is required to sustain 8500 maneuvers, the allowable stress for HPS is 
60 ksi. However, for gain stabilization, the stress can be increased to 70 ksi. This 
represents a 16% weight savings. 

Results for the 600”F analysis is shown in Figure 31. Similar observation can be 
made as in the previous case. Results again indicate a longer fatigue life for the 
conventional gain stabilization. Figure 32 combines results of both temperature settings. 
Note that elevated temperature exposure reduces the fatigue lives for both stabilization 
techniques. 

In connection with the reduced fatigue life for HPS, examination of the normal 
acceleration response can provide some insights. As shown in Figures 28 and 29, at the 
end of the maneuver, the normal acceleration for the HPS design displays a higher peak 
value at 15 seconds than that of the conventional gain stabilized design. This difference 
amounts to a 0.32g increase or a 16% change in peak value. In addition to the difference in 
peak acceleration value, the time history data for the Phase Ila HPS shows that a 2 Hz 
oscillation or residual response begins to develop at 3.5 seconds. Experience suggests that 
the peak normal acceleration value accounts for perhaps 75% of the contribution to the 
reduced fatigue life of a structure. The residual response accounts for roughly the 
remaining 25%. 

In the original analysis, the 15 second duration of the maneuver was chosen 
somewhat arbitrary under the assumption that this would not be critical to the comparative 
assessment of the two designs. However because of the long path lag to attitude changes, 
coupled with small differences in the closed loop pitch responses at maneuvering 
frequencies, the differences in peak normal acceleration result. In a refined analysis, the 
differences in steady state normal acceleration would be removed to make the maneuver 
comparable for both designs. Thus, the primary remaining difference would be the 2 Hz 
oscillation associated with the first structural mode which is not filtered in the HPS design. 
This should considerably reduce the disadvantage of the HPS design. 

Funher, the fatigue analysis is an assessment procedure, i.e., it is an analysis done 
after a control design cycle. Further development (discussed in Section VI-B) could 
re-optimizc the HPS design with explicit consideration of fatigue to further reduce or 
eliminate the disadvantage. 
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Figure 31. Ccmparisonof PrediaedFatig\ieLifeat600°F 
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Hgurc32. Coo^arison of Predicted Fatigue Life 




C . Human>centered Requirements 

Phenomena involving human pilots or passengers include "passive" considerations 
such as ride quality, but the focus here is on "active" considerations, i.e., manual control. 
Recent developments in hypersonic flying qualities (Ref. 8, 9, and 12) indicate that certain 
vehicle constraints unique to hypersonic aircraft (e.g., high minimum angle of attack limits) 
could necessitate complex and unconventional maneuvering. While such maneuvers might 
be routinely automated, provision for manual maneuvering at some level must be provided 
unless a radical departure from aircraft operational philosophy is accepted. That is, human 
pilots will be demanded by human passengers and adequate manual control provisions 
must be made for these human pilots. Manual control is, of course, a primary motivator 
for HPS in connection with effective time delay reduction, but any associated adverse 
effects on manual control must also be assessed. 

The primary concerns are vibration feedthrough and interaction with the pilot's 
neuromuscular system. Vibration feedthrough occurs when normal acceleration (or 
acceleration in other axes) at the pilot's station is transmitted through the pilot's 
neuromuscular system to the stick. Preliminary assessments can be made using data 
available from biodynamic experiments (Ref. 13). 

The phenomena can be treated with the model shown in Figure 33. The pilot/stick 
dynamics represent the composite of the seat, human pilot and stick mechanical 
characteristics that determine the transmission of structural vibration to the stick electrical 
transducer. A fly-by-wire system is, of course, assumed, but the stick may be either force- 
sensing or displacement sensing at this point. The command path gain/filter represents the 
dynamics between the stick electrical transducer and the pitch rate command signal. This 
includes stick filters generally used in fly-by-wine systems but not the mechanical stick feel 
elements which are included in the first block. The augmented vehicle dynamics are 
represented by the appropriate >^zp/qc transfer functions (Table 9). 

Figure 34 presents human pilot transmissibility data in the form of the frequency 
response of the stick force/pilot normal acceleration transfer function (in Newtons/g). The 
discrete data points (squares) were obtained from human vibration experiments (Ref. 13). 
The solid lines are Bode plots generated from a model fitted to the discrete data using the 
STI BIODYN software (Ref. 14). The major elements of the BIODYN model are depicted 
in Figure 35. The stick used to generate the Figure 34 data was a center stick with a very 
stiff spring restraint similar to that of some fly-by-wire controllers. While this controller 


62 




63 


Figure 33. Biodynamic Analysis of Manual Orntrol 




Figure 34. Frequency Response of Human Transmissibility 



Torso 


\X 1 


Rotations 


I<hn. 

Bhn 


•"'Cov^ 


Activo 

Neuromuscular 

System 


Control Slick 
Mass,Oomping 
and Compliance 




’WSg/^Upper Arm: Forearm: 


Interface 

Deflection 

. **1 . 


M,,I, 


Mq t Iq 



A 


V 


Control 

Deflection 

A 


5. Ly-I“* 

/b. 

Interface Stick 

Compliance Force 
and Transducer 

Damping 





Figure 35. Elements of BIODYN Biodynamic Model 


65 



may differ somewhat from what would be expected for future HSV applications, the data is 
considered to be adequate for an initial comparative biodynamic assessment of HPS and 
conventional designs. 

The presence of a detent/threshold around the zero position of a fly-by-wire stick 
will tend to reduce vibration transmission in cruise conditions. This factor will not be 
treated here, but the results obtained should be relevant to the maneuvering conditions of 
interest where the stick is generally out of the detent. 

1 . Pure Gain Approximation of "zp / qc 

The complete BIODYN model is relatively complex and high order and it is 
appropriate to make initial biodynamic assessments with much simplified models. As a 
first step, it can be noted from Figure 34 that the transmissibility magnitude is roughly 
constant over much of the frequency range of the aircraft structural modes. This would 
imply a pure gain first approximation except that there is a rapid phase angle decrease with 
frequency in this range. At this first level of approximation, some uncertainty must be 
accepted in phasing, i.e., in the sign of a pure gain approximation. Thus, both positive and 
negative gains were examined in the analyses. 

If a pure gain command path characteristic (no bandlimiting stick filter) is also 
assumed, then the Figure 33 system becomes a pure gain feedback loop closure around the 
augmented airframe tizp/qc transfer function. Figure 36 compares these pure gain closures 
as root locus surveys for the conventional gain stabilized design and the final Phase Ila 
HPS design. Loci are shown for both positive and negative gains to allow for phase 
uncertainty. Figure 37 shows the same closures (for positive gain only) in Bode root 
locus plots which give a clearer picture of gain sensitivity. It can be seen that for both 
designs there are several potential instabilities associated with several structural modes. 
The first instabilities (modes that go unstable at the lowest gain magnitude for both positive 
and negative gains) are indicated in each plot (hatched horizontal line segment). For either 
sign of the gain, a larger gain magnitude is required to destabilize the Phase Ha HPS design 
indicating an advantage over the conventional design. The minimum gains for instability 
are summarized in Table 1 1 . 
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Figure 36. Root Locus Survey of Pure Gain Biodynamic Coupling 
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Table 1 1. Summary of Biodynamic Coupling Comparison 




Minimum Gain (deg/sec/g) in ci^/n^p Eiement for Instabiiity 
Intogors in [] indicate structurai mode driven unstabie 


Model 

Gain Sign 

Desian 1 

Conventionai 

Reference 

Ph Ha HPS 



Pure Gain 

+ 

0.19 [2] 

0.22 [3] 


- 

-0.13 [3] 

•0.36 [4] 

Mechanicai 

+ 

0.25 [3] 

0.46 [2] 

Stick 

* 

-0.13 [2] 

•0.92 [2] 

Stick 

+ 

2.29 [2] 

1.19 [2] 

Fiiter 

- 

-3.43 [3] 

-0.81 [3] 


Estimated qg/^zp 


qc U qc 

n^ n^ f» 

Zp *p * 


0.728 °/sec/g for f^/n^ » 3ib/g 
0.078 Vsec/g for fg/n^ - 28ib/g 
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2. Reflned Transmissibility Model 

Examination of Figure 34 shows a high frequency magnitude rolloff beginning 
around 50-60 rad/sec. This implies that a pure gain approximation is overly conservative 
with regard to pilot-transmitted excitation of the higher structural modes. To address this 
and to better model the phase characteristics, a third order transmissibility model was used 
in a refined analysis. This more detailed model (Ref. 13) was obtained from a second set 
of data as shown in Figure 38. This figure shows an even lower breakpoint (35 rad/sec) 
for the transmissibility rolloff, which should make it less conservative than the Figure 34 
data. These dynamics are represented by the transfer function 

Js_ ^ 35^(-s + 80) ( 5 ) 

"zp “ (s 80) [s^ -H 0.03(35)s + 35^] 

Note that this transfer function has unity low frequency gain and thus it simply models the 
spectral shape of the transmissibility. 

Figure 39 shows the vibration feedthrough loop closures for this case. Again, the 
Phase Ila HPS design appears superior to the conventional design reflected in the higher 
stability margins. See Table 11. 

3. Fly-by-Wire System with Stick Filter 

As a final case, a controller configuration somewhat more representative of a modem fly- 
by-wire design was examined. This consisted of a relatively stiff stick (implying a high 
mechanical bandwidth) coupled with a stick filter (a command path electrical filter 
downstream of the stick transducer). This type of stick filter is largely intended to filter 
pilot remnant and thus tends to have a bandwidth well below the stick mechanical 
bandwidth. Therefore, the pilot/stick mechanical dynamics (transmissibility) was 
reasonably approximated by a (unity) pure gain. The command path dynamics created by 
the stick filter is modeled by 


9s. - 

fs " 


10 


(s + 10) 


( 6 ) 


Again, this transfer function has unity low frequency gain and represents the spectral shape 
of the command path. 
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Figure 38. Vibration Feedthrough Frequency Response for a Spring Stick 
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Figure 40 shows the normal acceleration feedback loop closures for this system as 
conventioniil nx)t loci and F'igure 41 shows the corresponding Bode root loci. Since the 
stick filter configuration has the lowest command path bandwidth of the three examined, 
the least excitation of structural modes would be expected and this is reflected in the 
Table 1 1 margins. Table 1 1 also shows that the gain stabilized reference design, with the 
stick filter, has higher margins than the Phase Ila design - a reversal from the situation in 
the two previous cases. With the stick filter, there is less uncertainty about phasing and die 
positive gain should be the physically meaningful case. On this basis, the conventional 
design has about twice the margin of the HPS design. 

Some insight into this reversal of phenomenon can be gained from comparing the 
Bode root locus plots of Figures 37 and 41. With the introduction of the stick filter, the 
higher structural modes become less critical and the margin at the first structural mode 
becomes comparable to that at the (critical) higher modes. In this situation, the gain 
stabilized design has an inherent advantage due to the "notching" of the first mode. 
However, the ultimate significance of this depends on a more fundamental question, i.e., 
the significance of the margins on an absolute basis. This is addressed next. 

4. Quantitative Assessment of Biodynamic Coupling 

The question remains as to whether there is an actual potential for a closed loop 
biodynamic instability for any of the SCAS design/conunand path configurations 
considered above. To assess this, some absolute values for gains in the Figure 33 system 
must be determined. This is summarized in Figure 42 and explained further below. 

The augmented vehicle dynamics have been given previously in Table 9. A representative 
magnitude for the pilot/stick dynamics in the structural mode frequency range can be 
obtained from Figure 34 as 30 dB in Newtons/g. This corresponds to 31.6 N/g or 
7.1 Ibs/g. 

The final element, the command gain, requires some development. For a steady 

pullup. 
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a) Gain Stabilized Reference Design 



b) Phase lla HPS Design 


Figure 40. Root Locus Survey of Biodynatnic Coupling with Fly-by-wiie 
Stick and 10 rad/sec Stick Filter 
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Figure 42. Quantificati(» of Bi<xlynaiDic Coupling 



For the Mach 6 flight condition, Vxo/g = 186.3 g/rad/sec or 3.25 g/deg/sec. 
Figure 43 compares this value with the low frequency asymptote of the **zp/qc frequency 
response and shows the actual system response to be close, but a little low. This is due 
largely to imperfect command following of the pitch loop which was not optimized. The 
difference, however, is insignificant for the purposes of this study. 

A point to be noted from the basic kinematics (Eqn. 7) is that hypersonic pitch rates 
must be much lower than fw conventional aircraft to avoid exceeding structural load liimts. 
However, it can be argued that the stick sensitivity, i.e. "stick force per g", will necessarily 
be much closer to conventional aircraft levels. For the purposes of this analysis, the 
MIL-STD-1797A Level 1 range of 3 to 28 Ibs/g will be used. According to the MIL-spec, 
the actual values depend on the parameter 

® 

As is now well recognized, HSV will have unconventionally low values of 1/Te2 
(Ref. 8) which in turn implies that the existing database on stick sensitivity may be 
inadequate or inappropriate for HSV design. However, the above noted stick sensitivity 
range should serve for an initial analysis here. 


Thus, given 


_ 9£ Dz „ (9) 

fs “ fs 9c g ) 


the low frequency gain of the command path element can be estimated as 



1.92x10“^ to 1.79x10'^ rad/sec/lb 
0.0100 to 0.103 7sec/lb 


( 10 ) 


It should be noted that the stick force per g transfer function fs/n^ is physically 
independent of the transmissibility transfer function fs/"zp. Also, under the short period 
approximation used in this work, the steady state is the frequency range relevant to 
maneuvering and thus to the stick force per g specification. 
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43. Qosed Loop Normal Acceleration Magnitude Plot at the Hlot's Station 
for Phase Ha HPS Design 



The low value in Eqn. 10 corresponds to the upper end (28 Ibs/g) of the stick 
force per g range. Conventional practice would call for higher levels at high speeds, but 
the low hypersonic nz/a clouds this issue. 

To complete this analysis, the low frequency gain of qc/”zp is calculated to 
compare with the values in Table 11. 

_ is_ 9c = I 0.728 7sec/g for fg/nz = 3 Ib/g „ 

"zp nzp fs X 0.078 7sec/g for fs/n^ = 28 Ib/g ’ ^ 

Comparison of Equation 1 1 with Table 1 1 (positive gains) indicates that only with 
the stick filter is there a stable margin, for either design, at low values of stick force per g. 
As noted above, the conventional design has about twice the margin of the HPS design, but 
whether either margin is adequate is somewhat uncertain. It should be noted that while 
positive margins imply stability, there is also an implication of some reduction in structural 
mode damping due to the pilot vibration transmission path, i.e., there is a potential for a 
special form of pilot induced oscillation. There is no existing speciHcation for this 
biodynamic margin. No comparable analysis for a conventional aircraft is readily available. 

As noted above, conventional practice would indicate higher stick force per g at 
higher speeds but this is associated with a corresponding increase in nz/a. The 
unconventionally low values of hypersonic nz/a create an uncertainty about the ideal HSV 
stick force per g and thus uncertainty about biodynamic coupling. There are other factors 
that may influence biodynamic coupling with stmcmral dynamics. In particular, command 
path filters (beyond simple stick filters) may be needed to "quicken" HSV path response 
which may in turn affect biodynamic coupling. Recent work (Ref. 9) indicates that special 
roll maneuvering requirements may lead to requirements for response types other than pitch 
rate command which further complicates the analysis problem. 

Given these uncertainties, it appears that these biodynamic issues should be 
considered in the stick path design for any HSV. Further investigation including 
specialized manned simulations are warranted particularly in connection with the 
development of HPS design criteria. The simulation problem for such an investigation has 
particular challenges which are addressed in the next section. 
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Section V 


APPROACHES TO VALH)ATION OF HPS DESIGN METHODS 


A concept such as HPS flight control design must ultimately be validated. In 
general, this implies manned simulation studies and flight tests. For HSV, this process is 
particularly challenging but necessary in proportion to the challenge. Detailed consideration 
of flight validation is premature in the sense that initial HSV flight programs, e.g. the 
NASP, will be research programs in themselves. However, it was appropriate to begin 
consideration, in Phase II, of manned simulation studies, regardless of whether the final 
Phase Ila HPS design is sufficiently developed to validate. This reflects the complex 
problems to be solved to do meaningful simulation. To a large extent, these problems are 
related to the very wide frequency range that must be simulated. However, the need for 
such wide frequency range simulation arises not only for validation, but for further 
understanding of the impact of residual response and development of HPS design 
requirements. Even beyond these particular motivations, there appears to be an inherent 
need for wide frequency range simulation in HSV development. For example, 
investigating such complexities as propulsion system/structural dynamic coupling. 

Of the issues addressed to date in this study, the questions of residual response arc 
most appropriately addressed with manned simulation. The issues of effective time delay, 
system stability and robustness can all be examined with manned simulations. These areas 
can also be addressed analytically much more readily than can residual response questions. 
And even in the area of residual response, it should be noted that only the human-centered 
issues are really appropriate to address with manned simulation. Problems such as 
structural fatigiw probably cannot be usefully studied directly with manned simulation. 

In the dichotomy of passive and active human-centered problems noted above, the 
passive problems of ride quality are certainly important for HSV and similar aircraft such as 
the high speed civil transport (HSCT). However, the data from manned simulation studies 
related to ride quality (Ref. 15 and 16) are much richer than that for active issues such as 
biodynamic vibration transmission. Because this problem has intrinsic complexities for 
manned simulation beyond those of ride quality studies, it was the focus of consideration in 
Phase II for developing manned simulation approaches. At this point, there is no real need 
to distinguish between the simulation requirements for HPS validation and the requirements 
for supporting HPS concept development which precedes validation. 
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The primary HPS-related issues for biodynamic simulation studies of HSV manual 
control inclu(te: 

• Overall assessment of the impact of residual response 

• Development of biodynamic coupling data 

Controller (stick) design, stick filter design and seat design 

• Impact of HSV maneuvers 

• Impact of high altitude disturbances 

• Impact of structural model uncertainty 

As noted above, the primary driver for simulator capability is the wide frequency 
range required for this application. The unconventionally low value of hypersonic 1/Te2 
will tend to drive the minimum frequency range associated with maneuvering below 
0.1 rad/sec. On the other end, simulation of the flexible modes could drive the simulator 
bandwidth above 10 Hz. Specialized simulators exist which can reach one end or the other 
of this frequency range, but generally not both ends. 

A brief survey of existing simulators was made in Phase II to identify possibilities 
for hypersonic research. The capabilities of one such facility, a relatively small moving 
base simulator with a custom built hydraulic hexapod base, are summarized in Figure 44. 
This figure shows both simulator limits and human subject tolerance limits based on 
established standards. It can be seen that simulation in this facility is limited by human 
tolerance rather than simulator capability. It is interesting to note that this simulator, which 
in many respects is among the most well suited for hypersonic research, was built by an 
agricultural implement manufacturer for tractor development! 

Existing moving base aircraft simulators designed for flying qualities and flight 
control research, such as the NASA LaRC VMS, can certainly accommodate the low 
frequency requirements. However, these simulators generally will not have the bandwidth 
desired for structural vibrations. This limitation is set primarily by high cab mass reflected 
in motion base lags and limiting although it might be possible to address the first structural 
mode (assumed around 2 Hz). Some sort of "shaker seat" might be obtainable to inject 
high frequency motion in series with the motion base to increase the effective bandwidth. 

There are specialized vibration simulators, such as the Air Force AMRL "six-mode" 
(Ref. 14), which can meet the bandwidth requirements for this application. However, 
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these facilities generally cannot achieve the low frequency requirements, i.e., they cannot 
simulate maneuvering. Further, there are typically limitations on mechanization of closed 
loop manual control experiments. For example, the AMRL six-mode could be run with the 
subject performing a tracking task, but for safety reasons the subject's input would not be 
fed back to the motion dynamics. During the Phase II work, discussions were held with 
AMRL personnel and the possibility of including the subject in the motion loop was 
discussed and it appears that there are possibilities for this. 

One possibility, suggested by the above simulator limitations, is the use of both 
moving base flying qualities simulators and vibration simulators in coordinated experiments 
that would cover the wide frequency range of interest. The challenge, of course, is in 
integrating the results from several experiments in different simulator types. One 
possibility involves the use of a biomechanical model in the form of appropriate software, 
such as the STI BIODYN package (Ref. 17), to interface data from several simulators. 
Among the capabilities of the BIODYN package is the ability to identify the model 
parametei^ (see Figure 33) using frequency domain identification procedures applied to the 
measured simulator data. 

To integrate data from several experiments in different simulator types, separate fits 
would be made with appropriate free parameters identified in BIODYN for each simulation. 
The aggregate of parameter values fitted would then constitute an overall model which in 
principle covers the frequency range of interest. This application of biodynamic modeling 
would require development and would be a significant effort in itself. However, it would 
provide a rational and systematic approach to the tricky problem of combining the results 
from distinctly different simulators. The end result would be a refined biodynamic model, 
tuned to the HSV situation, which could be used for development of physically-based 
design requirements or even used directly in the flight control design process. 
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Section VI 


CONCLUSIONS AND RECOMMENDATIONS 


A. Results and Conclusions 

The hybrid phase stabilization design procedures were refined in this contract. The 
final Phase Ila HPS design is a significant improvement over the Phase I HPS design. It 
has achieved greater reduction in effective time delay as well as residual response. The 
final Phase Ha HPS design also meets the MIL-F-87242 gain and phase margin (8 dB and 
60®, respectively) except for the third structural mode. This mode, however, is phase 
stable. 


The final Phase Ila HPS design has a 36 msec effective time delay which is 
41 msec less than Phase I HPS. This is also weU below the 218 msec value for the gain 
stabilized reference design which greatly exceeds the MIL-STD- 1797 A Level 1 maximum 
of 100 msec. 

The flexible HSV model was extended to provide normal acceleration at the pilot's 
station. Then, the Phase I residual response metric was revised by replacing the previous 
estimated model with the new model extension. Based on this revised metric, the residual 
response of the final Phase Ha HPS design was shown to be much improved over that of 
the Phase I HK design and equivalent to that of the gain stabilized reference (ksign. 

The sensitivity of HPS designs to structural uncertainties were assessed by 
evaluating the stability metric with variations in structural damping ratio, mode shape and 
frequency. For the purpose of comparison, the sensitivity of gain stabilized designs were 
also evaluated. Both techniques appear to be somewhat sensitive to variation in structural 
frequency. While gain stabilization tends to be more sensitive to lower fiequency variatiOT, 
phase stabilization can be mtwe sensitive to the higher fiequency variation. One technique 
does not have distinct advantages over the other. Both HPS and gain stabilization 
techniques have their own potential problems and unique robustness characteristics. 


The revised Phase II metric docs not appear adequate for design bas«i on the 
structural fatigue study and the analysis of biodynamic vibration feedthrough, but has 
potential for refinement based on such physical bases. 
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Assuming titanium 6-4 fatigue properties, the preliminary structural fatigue study 
predicts that while elevated temperature exposure reduces the fatigue lives for both the gain 
stabilized and HPS designs, gain stabilization exhibits a lower vehicle weight for a given 
number of maneuver cycles. However, re-optimization of the HPS design to explicitly 
accommodate fatigue could reduce or perhaps remove this shortcoming for HPS designs. 

Biodynamic analysis of vibration feedthrough was performed for the gain stabilized 
reference design and the final Phase Ila HPS design. The two designs are roughly 
comparable, but the gain stabilized design may have an advantage due to the notch filtering 
of the first structural mode. However, there is a potential for destabilization of structural 
modes by the pilot with either system. It was found that biodynamic analysis of manual 
control is quite sensitive to the assumptions about the controller (stick) and command path 
dynamics and the stick force per g level appropriate for HSV. 

Because its effective time delay is well below the MIL-spec requirement, the 
Phase Ila HPS design has apparent potential for further refinement to meet residual 
response requirements. The critical need to achieve this is refinement of the multiple 
aircraft and human-centered requirements in a form suitable for design. 

B . Recommendations for Future Research 

The final Phase Ila HPS design is a clear improvement over the earlier HPS 
designs, but still has some apparent shortcomings with respect to the conventional gain 
stabilized reference design. This is based primarily on the structural fatigue and 
biodynamic coupling assessments made after the last design iteration. However, this is not 
conclusive because these assessments are not validated design requirements. Furthermore, 
the conventional reference design fails to meet established design requirements (in 
particular, MEL-spec phase margin and time delay) so it cannot actually be considered 
superior to the HPS designs. 

The final Phase Ila HPS design satisfies the time delay requirement by a wide 
margin (36 msec which is much less than the 100 msec required) which suggests that this 
margin could be traded for improvements in other areas. This in turn suggests that the 
problem is reaching the stage where more formal design optimization methods may be 
appropriate to efficiently address increasing subtle tradeoffs among multiple requirements. 
To this point in the project, problem definition, basic control system architecture, and 
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physical interpretation have been emphasized to understand design requirements and their 
implications. Standard classical control system synthesis methods have been adequate and 
appropriate for this activity, and thus, synthesis methodology has not been emphasized. 
With considerable basic development work done, it now appears appropriate to begin 
considering how optimal synthesis methods could be exploited to carry the HPS 
developments forward. 

However, to lead to useful design procedures for hypersonic flight control, a 
"problem-driven", as opposed to "methods-driven", approach to optimal synthesis methods 
must be adopted. That is, the procedure must begin with the best possible development of 
the design requirements for the HSV control problem. A search of existing optimization 
methods can then be made for those with the best potential for this application. Next, 
refinement of selected synthesis methods to meet the specific requirements of this 
application and the HPS concept would be made. The design requirements of this 
particular application can be translated to meet the format of specific optimization 
algorithms as long as the intent of the original requirement is not compromised. The basic 
principal is that the method must accommodate the application and not the reverse. 

The basic elements of recommended future work include development of the 
following: requirements, basic architecture, design strategies, and synthesis methods. 

In terms of requirements, the premise is that formulation of a complete and 
appropriate set of requirements is the key step in the design problem. Without these, no 
optimal synthesis method can produce a good HSV design. TTie practical problem is that 
there are multiple requirements which may even conflict. Major requirements include: 
effective rigid body dynamics, effective time delay, control power, MIL-spec gain and 
phase margins, specialized robustness specification, human-centered residual response 
requirements, and aircraft-centered residual response requirements. 

Regardless of the requirements and synthesis procedures used, the basic system 
architecture must be maintained under the direct control of the designer, the methodology 
must react accordingly. This includes the imptwition, by the designer, of the standard filter 
forms such as notches and lags with constraints on order and location. For HPS design, 
this must also allow the blending of specified multiple sensor signals with constraints on 
location. Ultimately, the inclusion of estimators should be provided for. 
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For a given architecture, there are still degrees of design freedom that the designer 
must be able to constrain as part of a design strategy. For HPS designs, the major issue is 
deciding which modes will be primarily gain stabilized and which will be phase stabilized. 
A directly related strategy issue is the tradeoff between damping a structural mode and 
filtering it out. Part of the overall strategy used to date is to treat the structural modes in a 
way which does not significantly affect the effective rigid body mode. This has the 
significant advantage of allowing the rigid body compensation to be designed 
independently, but it also constrains the positioning of the first flexible mode zero in HPS 
designs. Furthermore, integration of the rigid and flexible control design might allow 
further improvement in the HPS model damping. The real issue here is to improve the 
designer's productivity by reducing his computational load without eliminating the added 
value of his design guidance. 

Efforts are being made to investigate applications of new methods such as 
H-infinity control to the control of flexible vehicles in ways which address the above 
issues. Progress has been made in guiding H-infinity designs procedures to produce 
standard lag and notch filters. More importantly, procedures for damping modes in 
H-infinity designs have been developed which are relevant to HPS applications. In 
addition to these H-infinity approaches, validated software, having the capability fw fixed 
form optimization, exists which could be used to c^timize HPS architectures. 
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Appendix A 


LOW FREQUENCY GAIN MARGIN 


The main focus of HPS design, once the superauginented pilch loop is established 
at the rigid body level, is on the high frequency (structural mode) region. However, as can 
be seen in Table 1, the low frequency gain margins of all four designs (at the phase angle 
crossover around 1 rad/sec) are below the MIL-F-87242 minimum. The low frequency 
phase crossover arises from the inherent conditional stabiUty of the superaugmented pitch 
loop design applied to a statically unstable airframe (Ref. 18). The final Hiase Ha HPS 
design comes closest to meeting the requirement and the gain stabilized reference design is 
the most deficient. These deficiencies were considered minor with reasonably 
straightfOTwaxd solutions, but no attempt was made to rectify the designs because this 
deficiency has very little significance with respect to the critical questions of HPS design. 

However, the possibilities for improving the low frequency gain margin were 
examined briefly at the outset of Phase H, based on the Phase I HPS design, to verify the 
above assumptions. Lag-lead filters were inserted in the pitch loop to boost the low 
frequency gain margin by extending the Ic/s asymptote below 2 rad/sec. This can be 
accomplished by placing the compensator lag at 1.5 rad/sec and the lead at 2 rad/sec. The 
effects of this compensator are shown in Figures A-1 through A-4 which should be 
compared to Figure 7. As can be seen in Figures A-1 (Bode root locus plot) and A-2 
(Bode magnitude and phase plots), the k/s asymptote has been extended below 2 rad/sec, 
thus increasing the low frequency gain margin above the required 6 dB. Figure A-3 
displays the root locus for this system. Notice that the compensation creates a potential for 
increasing the damping of the first flexible mode. Because of the gain limitation impUed by 
the desired rigid body bandwidth, this potential is not exploited in this design. 


The closed loop pitch rate response to a step pitch rate command is shown in 
Figure A-3. Comparison with the original Phase I HPS design (Rgure 13) shows that 
the modified system approaches the steady state value considerably faster. This difference 
results from increased closed loop frequency of the effective short period naode. That is, 
the presence of the lag-lead compensator has changed the bandwidth of the design. This is 
not surprising, but does imply that further design iteration would be required if the original 
dominant mode target was to be achieved. However, as noted above, the only purpose of 
this exercise was to demonstrate that the low frequency gain margin problem would 
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respond to standard treatment. This treatment was not applied to later designs since the 
is.suc is no! considered signifienni with respect to the objectives of the project. 
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Figure A- 1 . Bode Root Locus Plot of Phase I HPS Design with Low Frequency 

Lag-lead Compensator 
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Figure A-2. Bode Plot of Phase I HPS Design with Low Frequency 

Lag-lead Compensator 
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Figure A-3. Root Locus of Phase I HPS Design with Low Frequency 

Lag-lead Cwnpensator 
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Figure A-4. Closed Loop Pitch Rate Response at the Pilot's Station to a Unit Step 
Command for the Phase I HPS Design with Low 
Frequency Lag-lead Compensator 


A-5 



This page is intentionally left blank. 


A-6 


Appendix B 


SIGNIFICANCE OF NEGATIVE PHASE MARGINS FOR HPS DESIGNS 


Phase stabilized control system designs (or the phase stable regions of HPS 
designs) have an identifying characteristic in the phase curve of their frequency responses. 
This is the "sawtooth Bode" characteristic seen in the Figure 16 phase curve. One 
important consequence of this is that the open loop transfer function phase angle generally 
remains between -180® and +180“ in the phase stable region. This is distinct from the 
phase characteristic of gain stabilized designs (cf. Figure 5, noting that the phase curve is 
"folded" whenever it crosses the edge of the plot). In these designs, the phase angle 
typically crosses -180® just before the first structural mode and decreases continuously 
through the structural mode region. 

These differences in phase characteristics imply a need for some care in interpreting 
phase between conventional and HPS designs. Specifically, while a negative phase 

margin for a gain stabilized design generally implies an instabiliQ^, this is not the case for an 
HPS design. For example, the complete tabulation of phase margins for the Phase Ila 
HPS design (Table B-1) shows two negative phase margins associated with the first 0 dB 
gain crossings at the first two structural modes. HP S designs will tend to have more 0 dB 
crossings than gain stabilized designs, thus creating more phase margins. In particular, 
HPS designs create two phase margins for every sawtooth that peaks above 0 dB. This is 
illustrated in Figure 16. 


Table B-1. Phase Margin Tabulation for Phase Ha HPS Design 


Structural 
Mode No. 

1 

1 

1 

1 

2 

2 

3 

3 

Zero dB Gain 
Crossing (r/s) 

3.46 

12.6 

12.9 

36.3 

37.8 

52.3 

54.8 

I^ase Margin 
(deg) 

60.63 

-96.7 

96.9 

-166.7 

79.7 

142.6 

27.1 
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Figure B-1 shows a sketch of the Bode plot around a phase stable dipole. Note 
that phase angle peak implies net phase lead from higher frequency dynamics not shown. 
The sketch of the phase curve around the dipole shows that when the phase angle of Gql 
is, say, +90®, then the phase margin is 270° which is the same as -90° and phase stable. 
The important question is what does this really indicates about robustness? The phase 
margin should indicate if the system is robust with respect to phase uncertainties. For 
flight control design, particularly with the rigid body models ordinarily used, the concern is 
with unmodeled high frequency lags due to actuators, sensors and (gain stabilizing) 
structural filters (often not well defined in preliminary FCS design). The simplest 
approximation for these is a pure time delay. 

For the case of pure time delay, a -90° phase margin implies a much more robust 
design than one with a +90° phase margin. The question that remains is that of robustness 
with respect to unmodeled phase lead. For the conventional case of flight control design 
using rigid body models noted above, unmodeled dynamics would be very unlikely to 
produce phase lead. However, when working with more complete models including 
structural modes, and detailed sensor and actuator models, this assumption may not be so 
reliable; conceivably the unmodeled dynamics in more sophisticated models could result in 
net lead. 

Thus, the possibility that unmodeled lead could drive a positive phase "blip" 
unstable should be at least considered in an HPS design. Subtleties occur, e.g. a -90° 
phase margin is more robust to unmodeled lead than a -45° margin, but its the other way 
around for unmodeled lag. This raises new questions of what phase margins should be 
with possible unmodeled lead. As a minimum, the existing MIL-specs should be applied 
thoughtfully (as they always should be in any case!). 
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Figure B-1. Illustration of Negative Phase Margin for a Phase Stabilized Design 
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Appendix C 


ACCELEROMETER DYNAMICS IN HYPERSONIC FLIGHT 
OVER A SPHERICAL EARTH 


As noted in Ref. 8, "rotation of vertical" (R.O.V.) due to hypersonic flight over a 
curved earth (for nominally level flight) does not effect the force equations. Rather, the 
effects of "centrifugal relief enter the dynamics through the kinematic relations between 
angular velocity (with components P, Q and R referenced to inertial space) and the Euler 
angles (M^, 0 and d> referenced to local vertical). The dynamics of an accelerometer (on a 
rigid HSV) will be formally developed here to exhibit the R.O.V. effects. It will be shown 
that R.O.V. effects are of such low frequency that they do not affect HPS design. 

A. ACCELERATION EQUATIONS 

From Figure C-1, 


m(Afl«a + x) = -cx-kx + mg-tfo 


(C-1) 



Figure C- 1 . Freebody Diagram of Accelerometer Seismic Mass 
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where 


Aa = acceleration of case with respect to inertial space 

ea = unit vectOT in direction of accelerometer sensitive axis 

m = mass 

C = damping coefficient (c> o) 

k = spring constant (k > o) 

g = gravitational acceleration vector 

X = displacement of seismic mass 

Rearranging, 

-m(x + ^x + ^x) = m (Aa -g)- (C-2) 


At fiequencies well below the accelerometer natural frequency Vk/m, 

- ^ X » (Afl - g) • ea 


(C-3) 


The accelerometer signal Sa is 

Sa s - ^ X = fAfl - g; • Ca 


(C-4) 


This formulation is conventional and applies whether or not the earth is assumed 
flat. Further, g does not (and should not) contain any centrifugal relief corrections; these 
effects enter through the kinematics. Inverse square correction for altitude would be made 
directly to g, but this is small (a few percent) for flight within the atmosphere. 


If ^ and i are defined as the derivatives of $ with respect to inertial frame and 
at 

body frame, respectively, then from Figure C-3, the velocity equation of the accelerometer 
in body axis can be written as 




= Vc^+ p + iixp 


(C-5) 
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where p » 0. 


Since 



Inertial Reference Point 

Figure C-2. Accelerometer Location 


dQ ^ dp 

Aa “ yeg + ^xVeg + dt 

= Veg + ayVeg + (fl + i2xi2;xp + ax(axp) 

= Veg + iixVcg + ^xp + Qx(£2xp) 
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expressing Ag in the aircraft body frame yields 
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Neglecting vertical and lateral offsets in the accelerometer, y = z = o. 


Aa = 






U - RV + QW - (R^ + Q^)X 
V - RU - PW - Rx + QPx 
W - QU + PV - Qx + PRx J 


(C-8) 


Expression of the gravity vector in the body axes introduces the Euler angles 
referenced to the local vertical (not inertial space). 


g 


g 


-sin0 

sinO COS0 
COS^> COS0 


(C-9) 
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Assuming the vertical accelerometer sensitive axis is aligned with the z body axis, 
i*e«, 

1 

1 

II 


= W - QU + PV - Ox + PRx - g cosO COS0 

(C-10) 

B. PERTURBATION EQUATIONS 


The reference condition is constant altitude flight in the vertical plane over a 
spherical earth (great circle wbit) in a wings level attitude. Thus, W©, Pq. Ro. Qo. and 
are all zero. The perturbatitm normal accelerometer output, az is 

8z « w - QoU - qUo - qx + (g sinOo) 0 

(C-11) 

where lateral-directional perturbations are assumed zero (p = 
pitch rate is (Ref. 8) 

r = 0). The steady orbital 

Qo = 

(C-12) 

where 


Rg = Re + H 

(C-13) 

is the orbital radius and 


Uo = VToCOs0o 

(C-14) 

Further (Ref. 8), 


q s= e-AVi/Rg 

(C-15) 
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BuIVt - V , thus 


out. 


AVt = i (2 u + 2w) / V Uo^ + Wo^ = (u + w) / Vjo (C-16) 


Eliminating AVt, 




(C-17) 


Thus, 


% 


= w + ¥^u-0Uo-v^^ (u + w)-qx + (gsineo)0 


'^g 


VToRg 


= w - eUo - 0x + (g sin©o) 0 




t_ 


R.O.V. terms 


And for constant speed (u = 0, short period) approximation, 
az = w - 0 UO - 0x + (g sin0o) 0 


COS0O 


w + 


VtoR, 


w 






R.O.V. terms 


(C-18) 


x/VjoRg is much less than 1 and , and thus, the last R.O.V. term can be thrown 
Taking the LaPlace Transform gives the following expression. 


8z — 


f COS0( 


s - 


Rn 


w - 0 UO - 0x + (g sin0o) 0 
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There is only one R.O.V. term (cos0o/Rg) with an effective frequency of 
10'^ rad/sec associated with it. This term should have negligible effect at the structural 
frequencies and can be neglected. Thus, approximating in the frequency range of interest 
for HPS design, the accelerometer output is conventional which is consistent with the "flat 
earth" approximation. 

a* = w-eUo- ex + (gsin0o)0 
and 

q B» 0 
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